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Abstract 



in t 

We discuss the gravitational lensing of gravitational waves from merging 

On 



neutron star binaries, in the context of advanced LIGO type gravitational 
wave detectors. We consider properties of the expected observational data 
with cut on the signal-to- noise ratio p, i.e., p > pq. An advanced LIGO 
should see unlensed inspiral events with a redshift distribution with cut-off 



^ . at a redshift z max < 1 for /i < 0.8. Any inspiral events detected at z > z max 

should be lensed. We compute the expected total number of events which 
are present due to gravitational lensing and their redshift distribution for an 
advanced LIGO in a flat Universe. If the matter fraction in compact lenses is 
close to 10%, an advanced LIGO should see a few strongly lensed events per 
year with p > 5. 

PACS numbers: 98.62.Sb, 98.80.Es, 04.80.Nn 
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I. INTRODUCTION 



An advanced LIGO may observe gravitational waves produced as distant close neutron 
star binary pairs spiral into each other. During the last stage of inspiral the binary emits 
copious gravitational waves, with increasing frequency as the orbital period decreases, until 
finally the pair collides and coalesces. LIGO aims to detect the waves emitted during the 
last 15 minutes of inspiral when the frequency sweeps up from 10 Hz to approximately 10 3 
Hz 0. 

In this paper, we discuss the gravitational lensing of gravitational waves from merging 
neutron star binaries, in the context of advanced LIGO type gravitational wave detectors. 
Following Ref. f2|, we consider properties of the expected observational data with cut on 
the signal-to-noise ratio p, i.e., p > p . An advanced LIGO should see unlensed events 
with a redshift distribution with cut-off at a small redshift z max < 1 for h < 0.8 We 
argue below that there may be a significant number of inspiral events detected at z > z max 
which can be detected because they are magnified due to gravitational lensing. We compute 
the expected total number of events which are present due to gravitational lensing and their 
redshift distribution for an advanced LIGO, for plausible choices of cosmological parameters. 

The aforementioned frequency range over which LIGO can detect neutron star binary 
inspirals corresponds to a wavelength of gravitational waves from 3 x 10 4 to 10 2 km. This 
wavelength will be much smaller than the characteristic scales of gravitational fields the 
gravitational waves are likely to encounter as they pass between the neutron star binaries 
and the Earth. This means that one may treat the propagation of gravitational waves in the 
geometrical optics limit . In other words the gravitational lensing magnification will be the 
same as for optical light and one may use the standard formulae from optical gravitational 
lens theory. 
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II. OBSERVATION OF UNLENSED EVENTS 



Neutron star binary merger rate at redshift z per unit observer time interval per unit 
volume is n m = ho (1 4- z) 2 i](z), where no is the local neutron star binary merger rate per 
unit volume, (1 + z) 2 accounts for the shrinking of volumes with redshift (assuming constant 
comoving volume density of the merger rate) and time dilation, and rj(z) = (1 + z) 13 describes 
evolutionary effects. We use the "best guess" local rate density, h ~ (9.9 + 0.6 h 2 ) h x 
lO^Mpc^yr -1 ~ l(T 7 /i Mpc~ 3 yr -1 - §,§. 

In the last stage of a neutron star binary inspiral, gravitational radiation energy 
losses should lead to highly circular binary orbits. In the Newtonian/quadrupole ap- 
proximation, for a circular orbit, the rate at which the frequency of the gravitational 
waves sweeps up or "chirps", is determined solely by the binary's "intrinsic chirp mass", 
M.q = (MiM 2 ) 3 / 5 /(Mi + M 2 ) 1 / 5 , where Mi and M 2 are the two bodies' masses. For a binary 
inspiral source located at redshift z, the detectors measure M. = M.§{l + z), which is referred 
to as the observed chirp mass. For a given detector, the signal-to-noise ratio is 



di is our luminosity distance to the binary inspiral source. r and C(/max) depend only on 
the detector's noise power spectrum. The characteristic distance tq gives an overall sense 
of the depth to which the detector can "see". For advanced LIGO, r = 355 Mpc. C(/max) 
reflects the overlap of the signal power with the detector bandwidth (0 < ( < 1). For source 
redshift z, ( ~ 1 for 1 + z < 10 [2.8M /(Mi + M 2 )]. ( ~ 1 is a good approximation in the 
context of this paper. G is the angular orientation function, it arises from the dependence 
of p on the relative orientation of the source and the detector, < G < 4. Although G 
can not be measured, its probability distribution has been found numerically in Ref. [0, 
P (G, < G < 4) ~ 5G(4 - G) 3 /256, P©(6, G > 4) = 0. 

The luminosity distance di(z) = (1 + z) 2 dA(z), where <1a{z) is the angular diameter 
distance. In a flat Universe with a cosmological constant Q\ = 1 — Qq > 0, || dA.{z) = 




(1) 



cH -\l + z)- 1 Jo 2 dw [fi (1 + w) 3 + O a ] 



1/2 
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The number rate of binary inspiral events seen by a detector on Earth with signal-to- noise 
ratio p > Po per source redshift interval is 



dN NL (> po) 
dz 



4irn (cH x y 



<1>a(z) 
cHo 1 



[l + z)r}(z) 



Ce(x), 



(2) 



where C@(x) = f£° dQ Pq(Q) is the probability that a given detector detects a binary 
inspiral at redshift z with signal-to-noise ratio greater than p , it decreases with z and acts 
as a window function; Cq(x, < x < 4) = (1 + x) (4 — x) 4 /256, Ce(x, x > 4) = 0. x is the 
minimum angular orientation function 



x 



hA(r ,p ,M 
where we have defined parameter A as in Ref. | 



l + z) 



7/6 



A = 0.4733 



dA(z) 
cHn 1 



Mo 



(3) 
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(4) 



355 Mpc/ Vl-2M y 

Note that A absorbs the dependence on detector and source properties po, r , and M.$. 

Given the detector threshold in terms of the minimum signal-to- noise ratio po, the maxi- 
mum redshift of the source that the detector can "see", z max , is given by Eq.(|l]) with = 4. 
For advanced LIGO, ^ max < 1 for h < 0.8. [§,|J The redshift distribution given by Eq.@ 
terminates at z — z ma x- 



III. OBSERVATION OF LENSED EVENTS 

For a source at redshift z* > z max , we denote its signal-to- noise ratio without lensing by 
p(z*) [given by Eq. ([!])] . The source can not be detected in the absence of lensing. With 
lensing, its signal-to-noise ratio becomes 

p*(z*) = ^pp(z*), (5) 

where p is the magnification. The source can be detected if p*(z*) > p , with p denoting 
the detector threshold. 
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The probability of a source at redshift z being magnified by a factor greater than fi is 
P(> /i, z) = Tl(z) y 2 {p), for Tl(z) <C 1. r L (z) is the optical depth for gravitational lensing, 
and y 2 (p) — pT 2 for p ^> 1. For point mass lenses, we use 



1 



p > 



(6) 



The above equation leads to underestimation of y 2 (p) for p close to 1, which has negligible 
effect for our purpose. For a source at redshift z* > z max to be detected, we need \i > p , 
with 

2 



/'() 



Po 


2 












(7) 



where x(z*) is given by Eq.@. Note that /i depends on the angular orientation function 
0. The larger 0, the larger the signal-to-noise ratio without lensing [see Eq.(|l])], the smaller 
the magnification needed to reach the detector threshold po. Since = x is the minimum 
angular orientation function needed for a source to be seen without lensing [see Eq.(|3|)], 
only events with < x need be considered when we count the number of events which are 
present due to gravitational lensing. 

The number rate of binary inspiral events which can be seen due to gravitational lensing 
by a detector on Earth with signal-to-noise ratio p > po per source redshift interval is 



dN L (> Po) 



Ann (cHo 1 )' 



d A {z*) 



{1 + z*)t](z*)t l (z*) /■*(*•) 



d0P (0)y 2 M. (8) 



Note that m i n = 0, because we take the maximum magnification to be infinite, which is a 
reasonable approximation in the context of this paper. 

We consider two types of lensing: 1) "macro- lensing" from the large-scale gravitational 
field of galaxies, and 2) "micro-lensing" from the smaller scale gravitational field from com- 
pact objects such as stars. The optical depth from macrolensing is 
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T?(Z) 



F 
30 



1 + z) d A (z) 
cHo 1 



(9) 
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where F parametrizes the gravitational lensing effectiveness of galaxies [as singular isother- 
mal spheres]. Denoting the matter fraction in compact lenses as Q^, the optical depth of 
microlensing is [ |T0| . |TT| 1 



2 \(z) h J^(i + wf + n A 



where the affine distance (in units of cH x ) is A (2;) = J 2 dw (1 + w) 2 [fi (l + w) 3 + ^ 2 • 
For fi A = 0, rf (2) = § fi L [(y 5 / 2 + 1) lny/(y 5 / 2 - 1) - §], where y = 1 + 

IV. PREDICTIONS/SPECULATION 

To make more specific predictions we must choose parameters and to do this we are 
forced to speculate on the rate of inspiral events and the sensitivity of future gravitational 
wave detectors. Typically we expect neutron star binaries to have = 1.2 M while 
the advanced LIGO might have r = 355 Mpc, p > p = 5 then implies A = 0.7573. We 
parametrize evolutionary effects by r](z) = (1 + z) 13 , with f3 > and a redshift cut-off of 
z st0 p = 2.5. The typical lifetime of a neutron star binary is about 10 8 to 10 9 years ||; a 
significant fraction of neutron star binaries formed at z = 3 would have merged by z = 2.5. 



Since there seem to be a lot of star formation at z > 3 [12], -2 s top = 2.5 is probably reasonable. 

We find that the macrolensing rate is largest when there is a sizable cosmological con- 
stant, but is still negligible unless there is significant evolution. Even with extreme evolution, 



e.g. £7a = 0.8, fio = 0.2, F = 0.05 JTIJ and (3 = 3 macrolensing yields only 70% of the num- 
ber of lensed events of microlensing with more modest parameters: £7l = 0.07 and (3 = 0. 
Macrolensing of gravitational waves due to galaxies is negligible compared to a plausible 
microlensing rate. This is partly due to the fact that we have a good idea of the number 
and properties of galaxies while we are more free to speculate on the number of compact 
objects and partly due to the fact that point mass lenses are more effective gravitational 
lenses than galaxies. 

For the rest of the paper we restrict ourselves to microlensing, which gives a few strongly 
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lensed events per year without much evolution, for a currently acceptable value of Ql- We 
have considered two plausible cosmological models: (1) Q = 1 and = 0.1; (2) = 0.8, 
tt = 0.2, tt L = 0.07. 

We consider expected data with cut on the signal-to-noise ratio, p > 5. Fig.l shows the 
expected total number per year of events which are present due to gravitational lensing as 
function of h, for two cosmological models, with (3 = 0, 1. The solid lines are for Qq = 1 
and i\ = 0.1, and the dashed lines are for fl\ = 0.8, Qo = 0.2, and i\ = 0.07. Fig. 2 shows 
the corresponding expected total number per year of events which can be seen without 
gravitational lensing as function of h, with the same line types as Fig.l. The expected total 
numbers in both Fig.l and Fig. 2 increase with increasing j3, as expected. Fig. 3 shows the 
redshift distribution of expected events corresponding to Figs. 1-2 for h = 0.8. The dotted 
lines indicate the distribution of expected events which are present due to gravitational 
lensing. Note that gravitational lensing leads to tails at high redshift. For each cosmological 
model, the higher tail corresponds to /3 = 1. Note that in principle, the evolutionary index 
can be measured from the region of the redshift distribution dominated by events which 
can be seen without gravitational lensing. Note also that most of the events which are seen 
due to gravitational lensing lie beyond the cut-off redshift of the events which can be seen 
without gravitational lensing. 

We have used A4o = 1.2 M as the typical intrinsic chirp mass of neutron star binary 
inspirals. It is expected that AAq will fall in the narrow range of 1.12 — 1.26 M while 
an advanced LIGO can measure the observed chirp mass Ai = (1 + z)AA$ to an accuracy 
of better than 0.1% [0,0]. Thus the uncertainty in the redshift of a given event will be 



very small compared to the large range of z over which the events which are seen due to 
gravitational lensing are distributed [see Fig. 3]. Since the redshift distribution of observed 
events which can be seen without gravitational lensing should terminate at a relatively small 
redshift z max , an observation of an event with redshift significantly greater than z max is a 
strong evidence for gravitational lensing. One should be able to identify the events which 
are seen due to gravitational lensing! 



V. DISCUSSION 



While most neutron star binary inspiral events detected by an advanced LIGO will 
probably not be affected by gravitational lensing, there could be a detectable number of 
events which are significantly magnified via gravitational lensing by compact objects. These 
lensed events will be easily identifiable by their high observed chirp masses. For the no- 
evolution parameters used above one would expect around two events per year which are 
seen due to gravitational lensing. Even a modest evolution (/3 = 1) of the rate of inspirals 
can significantly increase the rate of events which are seen due to gravitational lensing, and 
one could imagine even stronger evolution. During the lifetime of a detector, say ten years, 
one might detect dozens of events which are seen due to gravitational lensing, from which one 
could estimate the amount of matter in compact lenses, Q^. The absence of such expected 
lensed events will place an interesting constraint on f^. 

If we can determine S\ m this way, we will have a much better handle on the nature 
of the dark matter in our Universe. Thus lensing adds utility to the observation of inspiral 
events, which has already been shown to provide a measure of the Hubble constant, the 
deceleration parameter, and the cosmological constant [ \A , . Gravitational lensing will 



also add additional noise to the determination of these cosmological parameters, although 
this noise is relatively small |T5|]. This is because, as we have seen, most inspiral events 



are little affected by gravitational lensing. Finally we note that our consideration of lensing 
for inspiral events is much the same as that which one uses when considering supernovae 



'standard candles" 16,17 
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Figure Captions 

Fig.l The total number per year of expected events which are present due to gravitational 
lensing as function of h, with (3 — 0,1. The solid lines are for Q = 1 and = 0.1, and the 
dashed lines are for Qa = 0.8, f2 = 0.2, and £\ = 0.07. 

Fig. 2 The total number per year of expected events which can be seen without gravita- 
tional lensing as function of h, with the same line types as Fig.l. 

Fig. 3 The redshift distribution of expected events corresponding to Fig.l for h = 0.8. The 
dotted lines indicate the distribution of expected events which are seen due to gravitational 
lensing. For each cosmological model, the higher tail corresponds to j3 — 1. 
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Figure 1 
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Figure 2 
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Figure 3 
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